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^ I Abstract: We investigate two models in non-commutative (NC) field theory by means 
of Monte Carlo simulations. Even if we start from the Euclidean lattice formulation, 
such simulations are only feasible after mapping the systems onto dimensionally reduced 
matrix models. Using this technique, we measure Wilson loops in 2d NC gauge theory 
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lO ! of rank 1. It turns out that they are non-perturbatively renormalizable, and the phase 
^ ' follows an Aharonov-Bohm effect if we identify 9 = 1/B. Next we study the 3d 

model with two NC coordinates, where we present new results for the correlators and the 
C^l , dispersion relation. We further reveal the explicit phase diagram. The ordered regime 
^ ' splits into a uniform and a striped phase, as it was qualitatively conjectured before. We 
also confirm the recent observation by Ambj0rn and Catterall that such stripes occur 
even in d = 2, although they imply the spontaneous breaking of translation symmetry. 
<D '. However, in d = 3 and d = 2 we observe only patterns of two stripes to be stable in the 



Oh! 



X 



range of parameters investigated. 

1 NC field theories as matrix models 

NC field theory is highly fashionable, in particular due to its potential role as a low energy 
effective description of string and M theory. Here we study Euclidean NC field theory as 
such; for a discussion of the Wick rotation and the related issue of unitarity, see e.g. Ref. 
|l|]. We assume two coordinates to obey the standard NC relation 

= iOef,^ . (1) 

This implies an uncertainty relation Ax^Axi^ = 0{6), and therefore UV and IR diver- 
gences are mixed. That property complicates perturbative renormalization drastically ^ . 
We can also observe non-perturbative manifestations of UV/IR mixing, see below. 

NC gauge theory can be formulated by a consistent use of the star product. Even 
for U{1) a YM term appears in the gauge action. 



S[A] = - I d^x Tr(F^, ^ F^,) , F^, = d^A, - d,A^ + ig{A^ ^A,~A,^ A^) . (2) 
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A lattice formulation can be written down, but its use in simulations is hardly possible: 
already the generation of star gauge invariant link variables appears as a fatal obstacle. 
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However, U{n) NC gauge theories can be mapped onto certain forms of the twisted 
Eguchi-Kawai model (TEK) ^ (which is defined on a single site), 

Stek[U] = -Nf3 Z,,Tt{U,UMIuI) , Z^, = Zl^ = e^-^/^ {fi<u, keN, L = N^'"), 

where f/^ are unitary N x N matrices and Z^^, is the twist. This mapping is based on 
"Morita equivalence", which is exact at ^ oo. At finite N there is an exact mapping 
onto a NC gauge theory on a lattice with spacing a in a finite volume 0, if we choose 
k = {N + l)/2, which implies 6 = No? /ti. After this mapping the system can indeed 
be simulated [^]. Thus we by-pass the problems with a direct lattice formulation. The 
limit N ^ oo corresponds to the simultaneous continuum and infinite volume limit (these 
limits are entangled due to UV/IR mixing). 



The action of the NC model reads 

"1 ^ , „ , ,0 A 



5[0] = y d^a; [-9^09^0+— 0^ + -0*0^0 . (3) 

Here the star product does not affect the bilinear terms, so the strength of the self-coupling 
A also determines the extent of NC effects. 

We consider d = 3 with a commutative time direction and two space coordinates 
obeying relation (|l|) Q Again we start from a lattice formulation — 4>{t,x) is defined 
on a T X A^^ lattice — and map it on a matrix model ||^ 0. Its action takes the form 

T 

= Tr E (r.'^(^) - kt))' + 1 {kt + 1) - + + ^0'(t) 
t=i ^ 

where 0(t) are Hermitian N x N matrices. The "twist eaters" provide the spatial 
shift, since they obey the Weyl-'t Hooft commutation relation F^Fi, = Zy^VyV^ (where Z 
is the same twist as in the TEK describing NC gauge theory). The final formulation has 
a remarkable similarity to the corresponding model on a "fuzzy sphere" 0]. 



2 Numeric results 



2.1 Gauge theory on a NC plane 



We simulated pure 2d NC U{1) gauge theory — for analytical work on this model, see 
Ref. 0. The large A^ limit at fixed (3 coincides with the planar limit of the commutative 



theory which was solved by Gross and Witten [|T3|. It exhibits an exact area law 
for the Wilson loop. From this connection we identify the "physical area" of one lattice 
plaquette as ~ However, in contrast to the planar limit we take the double scaling 
limit where A^ ^ oo, /5 ^ oo, so that 9 oc N/ (3 is kept constant. Thus we can extrapolate 
to a NC continuum theory in an infinite volume. 

The (analogue of the) square shaped Wilson loop in the TEK takes the form 

W^,{I xl) = ZZ TTiU^Uyuy) = W;,{I X I) , (4) 



^Including also the time in the NC relation seems more consistent from the gravity inspired motivation for NC geometry 
|lo[^, but it causes especially severe problems related to causality . 



corresponding matrix model formulation for 3d formions was suggested in Ref. [Q. 
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Figure 1: The polar coordinates of the complex Wilson loop W{I) plotted against the physical area A = 
a^I^. At small areas it follows the Gross-Witten area law (solid line). At larger areas the absolute value 
does not decay any more, but the phase grows linearly. 

and it corresponds to the Wilson loop in NC gauge theory. Note that the expectation 
value W{I) := j^{Wi2{I x J)) is complex. Fig. |l| shows its behavior in polar coordinates 
(for N/j3 = 32). The large N double scaling reveals that the Wilson loop is indeed 
non-perturbatively renormalizable. At small area (relative to 6) the result follows the 
Gross-Witten area law, hence in this regime the double scaling limit coincides with the 
planar limit of both, the commutative and the NC model. At larger areas, however, |iy(J)| 
does not decay any further, but the phase (p increases linearly in the area. Comparison 
of different values for N/ (3 shows in that regime the simple relation 

= A/9 = AB {A = {alf = physical area) (5) 

to hold to a very high precision. Here i? is a magnetic field across the plane, which we 
introduce formally as -B = 1/6', in agreement with the Seiberg-Witten map and with 
solid state applications |T^. This relation is illustrated in Fig. || (on the left). 

We also found large double scaling in some regime for the Wilson 2-point function 
and for the 2-point function of the Polyakov line (see Fig. ^ on the right), if a suitable 
(universal) wave function renormalization is applied 0. 
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Figure 2: On the left: The phase of the Wilson loop W{I) at N — 125 for various values of N/j3, 
corresponding to different values of 6. At large areas A > 0{6) the phase ip agrees with the formula 
Lp — A/ 9, which is shown by straight lines. On the right: The Polyakov 2-point function at N — 195, 
after performing a universal wave function renormalization, which provides a double scaling regime. 



2.2 The NC Xcp^ model in three dimensions 



The NC X(f)^ model has been studied analytically with various techniques. As a quali- 
tatively new feature (compared to the commutative world), it was conjectured that in 
part of the ordered regime (at strongly negative m^, which roughly corresponds to low 
temperature) stripe patterns may dominate over the Ising-type uniform ordering. Gubser 
and Sondhi performed a self-consistent Hartree-Fock type one- loop calculation , which 
would be exact in the 0{N) model at large N. They conjectured that it could also apply 
to the model, which suggests a striped phase in = 4 and d = 3, but not in d = 2 
(in the dimensions d < 4 they considered an effective action as suggested by Brazovskii 
|I7|). However, stripes did occur in 2d simulation results by Ambj0rn and Catterall ||T8| ; 
we add some remarks on that case in the Appendix. In d = 3 and 4, Gubser and Sondhi 
predicted a first order phase transition from the uniform to a striped phase if 6 increases 
in the ordered regime. At very strong 6 they expect more complicated patterns (like 
rectangles or rhombi) to become stable. 

Chen and Wu carried out a renormalization group analysis [|T^] in d = A — e. They 
predicted a striped phase for 9 > 12/y/e, which should therefore not occur in d = 4. 

We considered this model in d = 3 with two NC coordinates, and simulated the corre- 
sponding matrix model (described in Section 1) at T = = 15, 25, 35 and 45. Thus we 
obtained the explicit phase diagram shown in Fig. ^ It is based on the order parameter 



M{k) 



max 

NT \n\=k 



T 



Y.~^{t,p) 



t=i 



P=^n, ni,n2e{0,l,...,A^-l} . (6) 



{M{k)) detects a uniform order at /c = and a pattern of two stripes at A; = 1. Typical 
snapshots of ordered and disordered configurations at small and at large A are shown in 
Fig. Ij. We did not find more complicated patterns to be stable, but we consider it as an 
open question if this may set in at larger A^ or A (we recall that this corresponds to larger 
6) . The order-disorder transition seems to be of second order, whereas the uniform-striped 
transition appears to be first order — it takes the factor A^^ on the axes to stabilize it. 
In order to localize the transitions accurately we also looked for peaks in the connected 
2-point functions of M{k); some examples were given in Ref. 0. 
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Figure 3: The phase diagram of the 3d NC (j)^ theory. The connected points show the separation line 
between the disordered and the ordered regime, and the vertical lines mark the transition region between 
the uniformly ordered and the striped phase. 
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(b) Af2-m2 = -360 
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Figure 4: Snapshots of typical configurations (j)(t,x) at some time t, at N = 35, A''^A = 70 (a,b; disor- 
dered,uniform) and N'^X — 700 (c,d; disordered, striped). We show the {xi,X2) plane, where dotted and 
blank regions correspond to the different signs of (p. 



Here we illustrate the different phases further by showing the spatial correlator 



T 



C{x) 



t=i 



in Fig. ^. It does not decay exponentially in the disordered phase, see Fij 
Hence NC effects are visible also in the disordered phase. 

On the other hand, the temporal correlator in momentum space 



(7) 

on the left. 



G{r,p) 



TiV2 



t=i 



does follow an exponential decay (resp. a cosh behavior at finite T). Fig. ^ on the right 
is an example at p = 0. Varying p we can identify the dispersion relation E{p), which 
we show in Fig. ^ still in the disordered phase. It follows closely the linear behavior 
_ p'2 _|_ j^2^ (solid line, where Mgfj is a free parameter). The deviation at large 
momenta is a trivial lattice artifact. However, we also observe a systematic jump at 
p = 0, only in the case where we are close to the striped phase. This confirms the 
energy minimum at k = A^|p|/(27r) = 1, i.e. the trend towards stripes rather than uniform 
ordering. Such irregularities in the dispersion at significant 6 may be of importance; note 
that a 9 deformed dispersion relation has also been postulated for photons . 
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Figure 5: The spatial correlator C{x) defined in equation ^ at N'^X = 90 (a,b) and N'^X = 900 (c,d) 
at N = 45, against xt. The circles and squares correspond to two orthogonal directions; in particular in 
plot (d) the circles (squares) show the correlator parallel (vertical) to the stripes. 
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Figure 6: Correlation functions at N = 35 in the disordered phase (near the phase transition) at N^rri^ = 
—70, N^X = 490. We show the spatial correlator ^ on the left, and the temporal correlator in momentum 
space at p = on the right. 
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Figure 7: The dispersion relation E^{p^) in the disordered phase at N = 45. We are close to the uniform 
phase N'^rn^ = -22.5, N'^X = 90 (left) resp. close to the striped phase at N^m^ = -280, N^X = 900 
(right). The solid line represents the continuum behavior in a Lorentz invariant theory. 

3 Conclusions 

We presented non-perturbative results for two types of NC field theories. Such theories 
are non-local in the range of 0{\/6). We observe basic deviations from the commutative 
world rather at large scales though, which can be understood from the UV/IR mixing. 

Simulating the TEK we obtained results for 2d NC U{1) gauge theory. The scaling 
behavior that we observed is the first evidence for the non-perturbative renormalizability 
of a NC field theory. The Wilson loop follows an area law at small area. At large area 
it can be described by the Aharonov-Bohm effect if we identify B = 1/6, which appears 
here as a dynamical relation (on tree level it was used before in the literature). 

In the 3d model on the lattice we explored the phase diagram in the - A plane. 
In the ordered regime (at strongly negative m^) we found at small A a uniform order 
of the Ising type, as in the commutative case. At larger A — which amplifies the NC 
effects — we observed a pattern of two stripes to be stable for = 15 . . .45. It is an 
open question if multiple stripes etc. stabilize at larger values of A^ or A, as conjectured by 
Gubser and Sondhi. We observed the same behavior in d = 2, where we also verified up to 
A^ = 45 that only two stripes are ultimately stable (see Appendix). Still this implies the 
spontaneous breaking of translation symmetry, which is possible due to the non-locality 

In the ordered regime, the spatial correlations are dictated by the dominant pattern: 
uniform as in the commutative case, or striped with strong correlation in the direction 
of the stripes and anti-correlation vertical to them. ^ This agrees with the predictions 
in Ref. [IG]. In the disordered phase the spatial correlators deviate from the exponential 



^The stripes that we observed were always parallel to one of the axes. 



decay. However, the correlators in momentum space do decay exponentially in time for 
all momenta. 

That property allowed us to study the dispersion relation in the disordered phase: at 
small A also the dispersion behaves qualitatively as in the commutative model, but at 
large A there appears a jump at p = as a NC effect. Hence the energy minimum is at 
\p\ = 2'n'/N, in agreement with the trend towards two stripes. It remains to be clarified 
how this irregular dispersion behaves at large A^, and if it can be related to the expected 
non-linear photon dispersion at 6' > 0; the latter makes experimentalists hope for a precise 
measurement of 9 (for a discussion, see for instance Ref. [pT|] ). 
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A The NC A^^ model in d = 2 

The occurrence of stripes in the ground state implies the spontaneous breaking of trans- 
lation invariance. The case ci = 2 is particularly interesting in this respect. Gubser and 
Sondhi argued based on the Mermin- Wagner Theorem that stripes cannot be stable in 
d = 2. However, Ambj0rn and Catterall pointed out that this Theorem is in general 
not applicable to non-local theories, and in fact they did observe non-uniform patterns in 
their numeric results for the NC A0^ model in d = 2 [jl8[ (where the two coordinates obey 
the NC relation (P). 

We also performed simulations in d = 2, and we fully agree that also there stripe 
vacua exist. Starting from hot configurations and proceeding in Metropolis steps which 
propose a full Hermitian matrix to be added to 0(t) (with a coefficient tuned for a good 
acceptance rate), we find after about 500 steps a rich structure of patterns, where in 
general pi and p2 are both non-zero, and often larger than 2tt/N . However, if we continue 
to thermalize much longer we see that only the two-stripe pattern is ultimately stable, 
just as in (i = 3. Q This we tested again up to A^ = 45. The same can be seen if we switch 
to a more efficient algorithm, which updates the (conjugate pairs of) matrix elements one 
by one (this reduces the thermalization time by a factor oc 1/A^). Still the variety of 
meta-stable patterns may be of importance; we present some examples for them in Fig. 
B The Goldstone Theorem still holds for NC field theories |22], hence it is an interesting 



question to investigate the massless fiuctuations related to translation symmetry. 




10 15 20 25 30 



10 15 20 25 
Xi 



15 20 25 30 
Xi 



10 15 20 25 30 



Figure 8: Examples for meta-stable patterns in two dimensions at N = 35, A 
deeply in the striped phase. We show the same sort of maps as in Fig. H. 



0.43, = —2, which is 



^One might expect that the Theorem extends to theories, which are non-local in only a small range, which is invisible 
at very large scales. However, this argument is not adequate for NC models due to the UV/IR mixing. 

^However, we verified carefully, both in d = 2 and in d = 3, that those two stripes are really stable; one consistently runs 
into them starting from all sort of hot configurations, and they show absolutely no trend to disappear in very long histories. 
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